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Abstract 



o 
o 

■ In wireless networks, communication links may be subject to random fatal impacts; for example, sensor net- 



works under sudden power losses or cognitive radio networks with unpredictable primary user spectrum occupancy. 
Under such circumstances, it is critical to quantify how fast and reliably the information can be collected over 
attacked links. For a single point-to-point channel subject to a random attack, named as a dying channel, we model 
it as a block-fading (BF) channel with a finite and random delay constraint. First, we define the outage capacity 
as the performance measure, followed by studying the optimal coding length K such that the outage probability 



CS| ■ is minimized when uniform power allocation is assumed. For a given rate target and a coding length K, we then 

> 

ff^ I minimize the outage probability over the power allocation vector P^, and show that this optimization problem 

I ! 

. can be cast into a convex optimization problem under some conditions. The optimal solutions for several special 

(N 

. • ' cases are discussed. 

in 

. Furthermore, we extend the single point-to-point dying channel result to the parallel multi-channel case where 

' each sub-channel is a dying channel, and investigate the corresponding asymptotic behavior of the overall outage 

probability with two different attack models: the independent-attack case and the m-dependent-attack case. It can 

X 

' be shown that the overall outage probability diminishes to zero for both cases as the number of sub-channels 



increases if the rate per unit cost is less than a certain threshold. The outage exponents are also studied to reveal 
how fast the outage probability improves over the number of sub-channels. 
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I. Introduction 

Information-theoretic limits of fading channels have been thoroughly studied in the literature and to date 
many important results are known (see ^ and references therein). Generally speaking, if the transmission 
delay is not of concern, the classic Shannon capacity for a deterministic additive white Gaussian noise 
(AWGN) channel can be extended to the ergodic capacity for a fading AWGN channel, which is achievable 
by a random Gaussian codebook with infinite-length codewords spanning over many fading blocks such 
that the randomness induced by fading can be averaged out [|2l|[|3|l. With the transmitter and receiver 
channel state information (CSI) perfectly known, the adaptive power allocation serves as an effective 
method to increase the ergodic capacity. This allocation has the well-known "water-filling" structure O, 
where power is allocated over the channel state space. With such an allocation scheme, a user transmits at 
high power when the channel is good and at low or zero power when the channel is poor. When the CSI 
is only known at the receiver, the capacity is achievable with special "single-codebook, constant-power" 
schemes lIH. 

The validity of the ergodic capacity is based on the fundamental assumption that the delay limit is 
infinite. However, many wireless communication applications have certain delay constraints, which limit 
the practical codeword length to be finite. Thus, the ergodic capacity is no longer a meaningful performance 
measure. Such situations give rise to the notions of outage capacity, delay-limited capacity, and average 
capacity each of which provides a more meaningful performance measure than the ergodic capacity. 

In particular, there usually exists a capacity-versus-outage tradeoff for transmissions over fading channels 
with finite delay constraints |I71, where an outage event occurs when the "instantaneous" mutual information 
of the fading channel falls below the transmitted code rate, and a higher target rate results in a larger 
outage probability. The maximum transmit rate that can be reliably communicated under some prescribed 
transmit power budget and outage probability constraint is known as the outage capacity. In the extreme 
case of requiring zero outage probability, the outage capacity then becomes the zero-outage or delay- 
limited capacity [[8||. To study the delay-limited system, the authors in ^ adopt a i^-block block-fading 
(BF) AWGN channel model, where K indicates the constraint on transmission delay or the maximum 
codeword length in blocks. Such a channel model is briefly described as follows. Suppose a codeword 
is required to transmit within KB symbols, with the integer K being the number of blocks spanned by 
a codeword, which is also referred to as the interleaving depth (we call it coding length to emphasize 
how many blocks over which a codeword spans); it is also a measure of the overall transmission delay. 
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The parameter B is the number of channel uses in each block, which is called block length. A codeword 
of length KB is also referred to as a frame, where the fading gain within each block remains the same 
(over B symbols) and changes independently from block to block. The number of channel uses B in each 
block is assumed to be large enough for reliable communication, but still small compared to the channel 
coherence time. If the CSI for each J^-block transmission is known non-causally at the transmitter, transmit 
power control can significantly improve the outage capacity of the Ji -block BF channel flU. When the 
CSI can be only revealed to the transmitter in a causal manner, a dynamic programming algorithm is 
developed to achieve the outage capacity of the ii'-block BF channel in BH. 

In the above existing works, the delay limit is either infinite or finite but deterministic. However, 
there are indeed some practical scenarios where the delay constraint is both finite and random. For 
example, in a wireless sensor network operating in a hostile environment, sensors may die due to sudden 
physical attacks such as fire or power losses. Another example may be a cognitive radio network with 
opportunistic spectrum sharing between the secondary and primary users, where an active secondary link 
can be corrupted unpredictably when the channel is reoccupied by a primary transmission. How fast and 
reliably can a piece of information be transmitted over such a channel? This question motivates us to 
formally define the maximum achievable information rate over a channel with a random and finite delay 
constraint, named as a dying channel. This type of dying channels has never been thoroughly studied in 
the traditional information theory, and important theorems are missing to address the fundamental capacity 
limits. In this paper, we start investigating such channels by focusing on a point-to-point dying link and 
model it by a Ji -block BF channel subject to a fatal attack that may happen at a random moment within 
any of the K transmission blocks, or may not happen at all over K blocks. Note that the delay limit in the 
case of a dying channel is a random variable due to the random attack, instead of being deterministically 
equal to K as in a traditional delay-limited BF channel. Since the successfully transmitted number of 
blocks is random and up to K, a dying channel is delay-limited and hence non-ergodic in nature. Thus 
its information-theoretic limit can be measured by the outage capacity. It is well known that coding over 
only one block of a fading channel may lead to a poor performance due to the lack of diversity. However, 
when we code over multiple blocks to achieve more diversity in a dying channel, we must bear the larger 
possibility that the random attack happens in the middle of the transmission and renders the rest of the 
codeword useless. Therefore, it is neither wise to span a codeword over too many blocks nor just over 
one block. We need to consider the tradeoff between the potential diversity and the attack avoidance for 



the selection of the codeword length over such a dying channel. In other words, given a distribution of 
the random attack, we need to seek an optimal K that "matches" the number of surviving blocks in a 
probabilistic sense such that the achievable diversity is maximized and the outage probability is minimized. 

In a system with multiple parallel sub-channels (e.g., in a OFDM-based system), each sub-channel may 
be under a potential random attack. In such a scenario, we are interested in the overall system outage 
probability and how the outage probability behaves as the number of sub-channels increases. This leads 
us to examine the asymptotic outage behavior for the case of a parallel dying channel. We will consider 
two models of random attacks over the sub-channels: 1) the case of independent random attacks, where 
the attacks across the sub-channels are independently and identically distributed (i.i.d.) ; and 2) the case 
of m-dependent random attacks, where the attacks over m adjacent sub-channels are correlated and the 
attacks on sub-channels that are m-sub-channel away from each other are independent. 

In the following, we briefly summarize the main results in this paper: 

1) We introduce the notion of a dying channel and formally define its outage capacity. Suppose we 
code over K blocks, and the number of surviving blocks is random and up to K. An outage occurs 
if the total mutual information over the surviving blocks normalized by K is less than a predefined 
rate R. Correspondingly, the outage capacity is the largest rate that satisfies an outage probability 
requirement. 

2) We study the optimal coding length K that "matches" the attack time in a probabilistic sense 
such that the outage probability is minimized when uniform power allocation is assumed. We then 
investigate the optimal power allocation over these K blocks, where we obtain the general properties 
for the optimal power vector P^- We find that, for some cases, the optimization problem over 
can be cast into a convex problem. 

3) We further extend the single dying channel result to the parallel dying channel case where each 
sub-channel is an individual dying channel. In this case, we investigate the outage behavior with 
two different random attack models: the independent-attack case and the m-dependent-attack case. 
Specifically, we characterize the asymptotic behavior of the outage probabilities for the above two 
cases with a given target rate. By the central limit theorems for independent and m-dependent 
sequences, we show that the outage probability diminishes to zero for both cases as the number of 
sub-channels increases if the target rate per unit cost is below a threshold. The outage exponents 
for both cases are studied to reveal how fast the outage probability improves. 
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The rest of this paper is organized as follows. Section |ll] presents the system model for a single dying 
channel, as well as the definition of the corresponding outage capacity. In Section Hill we study the 
optimal coding length by considering uniform power allocation and derive the lower and upper bounds of 
the outage probability. Moreover, we obtain the closed-form expression of outage probability for the high 
signal-to-noise ratio (SNR) Rayleigh fading case. In Section |IVl we optimize over the power vector to 
minimize the outage probability. In Section |Vl we extend the single dying channel model to the parallel 
dying channel case. In particular, we examine the corresponding asymptotic outage probability with two 
setups: the independent-attack case and the m-dependent-attack case, in Sections IVTl and Ivnl respectively. 
In Section rvnn the outage exponents for both cases are examined to reveal how fast the outage probability 
improves. Section |IX] concludes the paper. 

Notation: we define the notations used throughout this paper as follows. 

• M indicates the set of real numbers, R_|_ is the set of nonnegative real numbers, and is the set of 
A^-dimensional nonnegative real vectors. 

• The error function: erf(a;) = e~^^ dt. 

• The normalized cumulative normal distribution function: = -k= C e~^dt. 

• The Q-function: Q{x) = e~~dt. 

• log(x) is the natural logarithm. 

• [■] is the ceiling operator and [-J is the flooring operator. 

II. Outage Capacity Definition of a Single Dying Channel 

We consider a point-to-point delay-limited fading channel subject to a random fatal attack, while the 
exact timing of the attack is unknown to neither the transmitter nor the receiver. Only the distribution of 
the random attack time is known to both the transmitter and the receiver. We further assume that there 
is no channel state information at the transmitter (CSIT) while there is perfect channel state information 
at the receiver (CSIR). The transmitter transmits a codeword over K blocks within the delay constraint; 
and when the fatal attack occurs, the communication link is cut off immediately with the current and rest 
of the blocks lost. We build our model of such a dying link based on the Ji -block BF-AWGN channel 
[[51, which is described as follows. 

Let X, y, and z be vectors in M^^^ representing the channel input, output, and noise sequences, 
respectively, where z is the Gaussian random vector with zero mean and covariance matrix a'^l^B- 
Rearrange the components of x, y, and z as K x B matrices, denoted as X, Y, and Z, respectively 
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(Each row is associated with B symbols from a particular block.). A codeword with length KB spans K 
blocks and the input-output relation over the channel can be written as follows: 



Y = AX + Z, 

where A = diag(|/ii|, ■ ■ ■ ,\hK\) is a K x K matrix with the diagonal elements being the fading amplitudes. 
Let Xj be the i-th column of X for z G {1, ■ ■ ■ , B}. Similarly, let Yj and Zj be the i-th columns of Y 
and Z, respectively. These are related as: 

Y, = AX, + Z,, ^G{l,---,5}, (1) 

which implies that the input symbols on the same row of X experience the same fading gain, i.e., they 
are transmitted over the same block. Since Z/s are i.i.d random vectors, we can view this channel as 
K independent parallel channels with each channel corresponding to a block. Hence, KB uses of the 
original channel corresponds to B uses of the K parallel channels in ©• The parallel channels over 
which a codeword is transmitted are determined by the channel state hi,h2, - ■ ■ ,hK, which can also be 



viewed as a composite channel lISl lfTOllH that consists of a family of channels {T{9),6 G 9} indexed 
by a particular set of 0. For a block fading channel with delay constraint K, it can be modeled as a 
composite channel {T(9) : 9 G ©a } as follows: Let Qk C be the set of all length- A' sequences of 
channel gains Ok = {hi,h2, - ■ ■ , /ix}, which occurs with probability tt^ under the joint distribution of 
{Hi, H2, ■ ■ ■ , Hk}- For each 0^ = {hi, /i2, ■ ' ' > ^x} ^ ©x^ we associate a channel T{0k), where T{Ok) 
consists of K parallel Gaussian channels. Let cxk = {ai,a2, ■ ■ ■ ax} be the fading power gain vector, 
i.e., ai = \hi\'^,i = 1, . . . , K, and = {Pi, P2, ■ ■ ■ , Pr} be the transmit power allocation vector. For a 
given set of Ok and Pk, the maximum average mutual information rate over channel T{Ok) is [HI: 

1 ^ 

Cbf{Ok, Pk, K) = —Y, log(l + (2) 

i=l 

where we assume a unit noise variance throughout this paper. 

In our model of the dying channel, the delay constraint is random rather than deterministically equal 
to K due to the fact that a random attack may happen within any block out of the K blocks or may not 
happen at all within the K blocks. If the fatal attack happens during the transmission, the current block 
and the blocks after the attack moment will be discarded. An outage occurs whenever the total mutual 

'a composite channel is a compound channel with prior probabilities. 
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information of the surviving blocks normalized by K is less than the transmitted code rate. Therefore, 
the dying channel is non-ergodic and an appropriately defined outage capacity serves as the reasonable 
performance measure. 

Let T be the random attack time that is normalized by the block length. As we know from the results 
of parallel Gaussian channels [[TTIl . with random coding schemes, we can decode the codeword even if 
the attack happens within the K blocks as long as the average mutual information of surviving blocks is 
greater than the code rate R of the transmission, i.e., if we have 

1 ^ 

-^log(l + «,P,)>i?, 
1=1 

then the codeword is decodable, where the random integer L = mm{K, [T\) with [-J being the flooring 
operator. 

Hence the outage capacity of a dying channel can be formally defined as follows: 
Definition 1: The outage capacity of a fT-block BF-AWGN dying channel with an average transmit 
power constraint P and a required outage probability 7] is expressed as 

CoutiP,v) = max sup <R: 

1 ^ 1 

^'^K ^ + "'^'^ <R}<v\- (3) 
i=i ^ 

Note that the outage probability above is defined over the distributions of the a/s and T, where we assume 

that the a^'s and T are independent of each other and the transmitter does not know the values of the a/s 

and T a priori, but knows their distributions. As we see from ([3]), there are two sets of variables to be 

optimized: One is the number of coding blocks K, and the other is the power allocation vector Pk- From 

the perspective of optimal transmission schemes, the outage capacity maximization problem is equivalent 

to the outage probability minimization problem dH . In the next section, we first study the optimal coding 

length K to "match" the attack time in a probabilistic sense such that the outage probability is minimized. 

in. Optimal Coding Length with Uniform Power Allocation 

As discussed before, we can optimize over the coding length K and the power vector to achieve the 
maximum outage capacity (or equivalently the minimum outage probability). If a uniform power allocation 
strategy is adopted, the only thing left for optimization is the coding length K. On one hand, we can have 
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a larger L = mm{K, [TJ) by increasing K, meaning that we potentially have higher diversity to achieve 
a lower outage probability. On the other hand, a larger K incurs a higher percentage of blocks being lost 
after the attack such that the average achievable mutual information per block is lower, and hence results 
in a larger outage probability. Since the random attack determines the number of surviving blocks and 
K determines the average base, we are interested in finding a proper value of K to "match" the random 
attack property in the sense that the outage probability is minimized. 

With uniform power allocation, according to the law of total probability, the outage probability can be 
rewritten as a summation of the probabilities conditioned on different numbers of surviving blocks, i.e.: 



where Aj- = {T,Li'^og{l + a iP) < KR} for j = !,■■■ ,K,Wi = Pr(z < T < z + 1) for i = 0, ■ ■ ■ , /sT- 1, 
and w*j^ = Pr(T > K). Given the distributions of and T, in general, there are no tractable closed-form 
expressions for Pr{Aj}'s. Alternatively, we could first seek the bounds of the outage probability and then 
study more exact forms for some special cases where we show how to find the optimal K. 

A. Outage Probability Lower Bound 

Notice that the following relationship holds: 




Wo + Pt{Ai}wi + Pt{A2}w2 + ■■■ 



+ Pt{Ak-i}wk-i + Pt{Ak}wI 



(4) 






> 



Pr <^ log (1 + aiP) < 



Since the fading gains a/s of different blocks are i.i.d, we have 





(5) 



where F{x) is the cumulative distribution function (CDF) of the random variable a,. 



Therefore, with the relationship in ([5]), we have a lower bound for the outage probability in (H]) as 



e — 1 ^ ' 



+ ^K- (6) 



B. Outage Probability Upper Bound 

On the other hand, there exists a simple upper bound for the outage probability: 



+ "^^) < ^ 1 < n + "*^) <KR}, j = l,---,K, 

I 1=1 ) i=l 

hence yielding 

i 

Therefore, an upper bound for the outage probability in ([4]) is given as 



'^^\^^Y.^og{l + a,P)<R^ 



+ [^[-^^]} <• (7) 

C. High SNR Rayleigh Fading Case 

From the previous discussion, we know how to bound the outage probability in terms of K with the 
general SNR values. However, there usually exists a significant gap between the lower and upper bounds. 
Fortunately, with appropriate approximations in the high SNR regime P for Rayleigh fading, we can obtain 
a tractable expression for the outage probability and hence further derive a closed-form solution for the 
optimal K. 

For our i^-block fading channel model with high SNR values, outage typically occurs when each sub- 
channel cannot support an evenly-divided rate budget (see Exercise 5.18 in ffT3l ). Thus, conditioned on 

^Here by high SNR, we mean that P is large, i.e., P S> 1. 
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the attack time T, the outage probability can be written as : 



Pout\T = Pr|-^^log(l + aiP) < i?| 



K 



Pr|log(l + a,P) < -i?|J . (8) 

For Rayleigh fading, we have Pr{ai < 1/x) ^ 1/x when x is large. Thus, when SNR is high, we can 
simplify ^ as 

Pout\T ~ (9) 

With the conditional outage probability given by dH), the overall outage probability is 

K 

Pout{K) = Wq + ^Pout\T ■ p{L = i) 



i=l 

= ^0 + 5^ ^^^i + (10) 

i=l 

Let G{t) be the CDF of the attack time, which is assumed to be exponentially distributed with parameter 

A. Let w*K = l- G{K), Wi = G{i + 1) - G{i) = - e"^) = p'c (for Vi < K) with c = 1 - e"^ and 

/3 = e~^. We can rewrite (flOl) as 

pUK) = e'"'Y,''-^ + ''-^\l-G{K)]+w„ 

= «'^«'^^ + ^^ + »o. (11) 
For high SNR, with < /? < 1, | is small. Hence, '^'^;'}^IP'' ^ when K > 2, and O can be 



^KR ^KR 



approximated to: 



PoutiK) ^ ee^^ + pj^J^ + wo, (12) 



where ,^ = (1 — e~^) -^J^ip - In order to obtain the optimal K by minimizing poutiK), we first treat (fT2l) 
as a continuous function of K, although K is an integer. 

Let us first consider the convexity of (fT2l) over a real-valued K. By taking the second-order derivative 
of (fT2l) over A', we have the following: 



d'poutjK) ^ [\ + log P-R] 



iR'e--+ ^-\:x_R,r ■ (13) 
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Since we liave A > and 1 — /5/P>0in tlie liigli SNR regime, it holds that ^ > 0. Therefore, (fT3l) is 
non-negative in the high SNR regime, which means that (fT2l) is convex over real-valued K. 

Given the convexity of (fT2)) . the optimal K can be derived by setting its first-order derivative to zero 
and finding the root. Consequently, the optimal solution K* is obtained as follows: 



K* = log 



X + logP-R 



A + 1okP 



Obviously, K* is unique given a set of ^, P, R, and A. Since a feasible K for the original problem should 
be an integer, we need to choose the optimal integer solution from [K*\ and [A'*], whichever gives a 
smaller value of (fT2l). 
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Fig. 1. Outage probability vs. coding length K, P=20dB. 

In Fig. \T\ and Fig. [2l we plot the coding length K versus the outage probability. We assume that 
the fading is Rayleigh, the random attack time is exponentially distributed with parameter 1/A = 10 
(normalized by the transmission block length), the target rate is i? = 1 nats/s/Hz, and the transmit power 
is set as 20 dB and 30 dB, respectively. As shown in Fig. [H the dashed curve and the solid curve are 
the lower and upper bounds given by ([6]) and dV]), respectively. The circles are obtained by using (fT2l) . 
As can be seen, firstly, the high-SNR approximation in (fT2l) is quite accurate. The circles are located 
between the upper and the lower bounds except for K = 1. This is due to the fact that when K = 1, 
^^^i~p/P ~ i-^P/P '^^^^ \:^o\dL . Secondly, we see that there exists a minimum outage probability 
over K as shown in Fig. |2l At last, comparing Fig. [Hand Fig. [2l we see that the upper and lower bounds 
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Fig. 2. Outage probability vs. coding lengtii K, P=30dB. 

get closer as the SNR increases with the values from (fT2l) are in between; hence the approximation in 
(fT2l) becomes more accurate. 

D. Low SNR Regime with Arbitrary Fading 

When SNR is low, we have log(l + ajP) ~ a^P. Thus, when we span a codeword over K blocks, the 
outage probability conditioned on T is given as 



Comparing (fT4l) and (fT5l) . we see that the outage performances of these two schemes are the same in the 
low-SNR regime. This is due to fact that in low SNR regime it is SNR-limited rather than degree-of- 
freedom-limited such that coding over different blocks does not help with decreasing the outage probability. 
Hence, repetition transmission is approximately optimal for a dying channel in the low SNR regime. 




(14) 



When using a repetition transmission (over blocks), the outage probability is given as 




(15) 
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IV. Joint Optimization over Coding Length and Power Allocation 

In the previous section, we investigated the optimal coding length K that minimizes the outage 
probability by assuming uniform power allocation. We now consider optimizing over both the coding 
length K and the power vector to minimize the outage probability. We note that optimizing over K 
is in general a 1-D search over integers, which is not complex. Since the main complexity of solving ([3]) 
lies in the optimization over P^, we first focus on the outage probability minimization problem over P^ 
for a given fixed K, which is expressed as: 



mm 

Pk 



K 

i ■ 

S.L 

K 

i=l 



^x:p.<p. (16) 



After obtaining the optimal outage probabilities conditioned on a range of K values, we choose the 
minimum one as the global optimal value. 

A. Properties of Optimal Power Allocation 

We start solving the above optimization problem by investigating the general properties of the optimal 
power allocation over a dying channel for a given K. 
Let Ej be the event that 

j 



K 

i=l 



It is obvious that the events are decreasing events, which means Ei ^ E2 ^ ■ ■ ■ ^ Ek- With the 
law of total probability, we can expand the outage probability in the objective of (fT6l) as follows. 



t{K) = Pr|-i^log(l + a,P,) <i?| 



Pout[ 

= wo + Fr{Ei}wi + Fr{E2}w2 + --- 

+ Fi{Ek^i}wk^i + Ft{Ek}w*k, (17) 

where w/s are defined in Section |llll With the above result, we then discuss the optimal power allocation 
for a dying channel under different conditions. 
1) Optimal Power Allocation over i.i.d. Fading: 
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Theorem 1: When fading gains over blocks are i.i.d., the optimal power allocation profile is non- 
increasing. 

Proof: The proof is provided in Appendi^iTAl ■ 
This is a general result regardless of the specific distributions of fading gains. That is, the optimal power 
vector lies in a convex cone !>+ = {Pr ^ ^+ ■ Pi > P2 > ■ ■ ■ > Pr}, matter what distribution the 
fading gain follows, as long as the i.i.d. assumption holds. 

2) Optimal Power Allocation over Identical Fading Gains: Now we consider the case where the fading 
gains over all the blocks are the same, while they are still random. This represents the case where fading 
gains are highly correlated in time. 

Theorem 2: When the fading gains a^'s are the same, the optimal coding length is K = 1 with Pi = P. 
Proof: The proof is provided in AppendixEB ■ 
This assertion implies that the optimal transmission scheme for a highly correlated dying channel is to 
simply transmit independent blocks instead of jointly-coded blocks. 

B. Power Allocations for Some Special Cases 

When the fading gain falls into some special distributions, we can further convert the corresponding 
optimization problem into convex ones and derive the optimal power vector efficiently. 

1) Optimal Power Allocation over i.i.d. Rayleigh Fading in High SNR Regime: Given ([8]) and condi- 
tioned on the attack time T, the conditional outage probability can be written as: 

L 



Pout\T = Pr < XI + aiPi) < KR 

I i=l 

^ f[Pr|log(l + a,P,) < (18) 
i=i ^ ^ 

For Rayleigh fading, we have Pr(«i < l/x) ^ 1/x when x is large. Thus, when SNR is high, we can 

simplify (fTSl) as 

PoutlT ^ J . (19) 

1 1.1=1 » 

The outage probability with Rayleigh fading in high SNR is approximated as below by substituting 
([H) into ([HI): 

PoutiK) ^wo + + ^ ' W2 + ■■■+ ^ J, ' <■ (20) 

^1 ^1^2 11^=1 Pi 
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Denoting Cj = Wi{e^^^^ — 1)*, we further simplify (|20l ) as 



Pout{K) ^wo + — + 



Ck 



Since the optimal power vector lies in a convex cone as shown in Theorem [H the problem can be 
formulated as a convex optimization problem (refer to AppendixEO for the convexity proof): 



Ci C2 Ck 

mm W0 + — + + ■ ■ ■ + 



.t. ^P, <KP, 



(21) 



i=l 



where = {P E : Pi > P2 > ■ • ■ > Pr > 0} is a convex cone. Thus, the optimal power vector 
can be efficiently solved with standard convex optimization algorithms such as the interior point method 

m. 

The simulation results are shown in Fig. |3l where we set the simulation parameters as: R = 0.5 
nats/s/Hz, 1/A = 5 for the exponential random attack, and average power P = 10 dB. As we can see, 
the power vector derived by solving problem (|2TI) achieves better performance in terms of the outage 
probability than the uniform power allocation case. 

1/X^5. P^IOdB 



- optimal power allocation 

- uniform power allocation 




2 2.5 3 

Number of blocks K 



Fig. 3. Outage probability with non-uniform and uniform power allocation. 



2) Optimal Power Allocation over i.i.d. Log-normal Fading: When the fading gain has a log-normal 
distribution, we can also approximate the problem as a convex one by minimizing the upper bound of the 
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objective function. Since we have 

L 



^log(a,P,)<I]log(l + a,P,; 



i=l 1=1 

the outage probability is upper-bounded as follows: 

L 



FiiJ^hgil + aiP,) <KR} 

i=l 
L 

< Pr{^log(a,P,) < A'i?} 



i=l 

K n n 

Wo + Yl ai<KR-Y^ log Pi}wn. 

n=l i=l i=l 



Thus, the optimization problem of (1161) is translated into the following problem, where we essentially 
minimize the upper bound: 

K n n 



min Wo + ^ P^^l^^ log < ^ log Pi}wn 

n=l i=l i=l 

K 

s.t. ^Pi< KP. (22) 



Let the ctj's be independent and log-normal random variables, i.e., logaj ~ A/'(0, 1), Vi. Since the sum 
of n standard normal random variables is a Gaussian random variable with zero mean and variance n, 
we have 

Pr{5^1oga, <x} = -{1 +erf(^)}, (23) 



1=1 * 
where erf(a;) is the error function. Substituting (l23l) into (|22|) yields the new objective function pout- 

K 

n=l * 

In general, (|24|) is not a convex function. However, under some special circumstances as described 
in Appendi}fDl the problem in (l22l) with the objective replaced by (|24|) can be rewritten as a convex 



Pout = wo + f2l{l+ erf( ^^ ^^^"^^^ )}^n. (24) 
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problem, which is given as following: 



mm Wo 

Pk 

K 



+ erf( -L^ )}Wr, 



1=1 

KR-logPi < 

/sTi?- log Pi - log Pa < 



K 



KR-^\ogPi < 0. 



(25) 



Therefore, efficient algorithms can be applied to solve the above problem. 

Numerical results are provided as follows. Assume that the outage probability target is set as ?7 = 0.3, 
the attack time is an exponential random variable with parameter 1/A = 4, and the fading gains are 
standard log-normal random variables. As we see from Fig. |4l the optimal power allocation leads to a 
significantly larger outage capacity over the uniform power allocation case. Moreover, as K increases, 
the outage capacity with the optimal power allocation may even increase to a maximum value while 
the outage capacity with uniform power allocation monotonically decreases. This suggests that, with the 
potential of a random attack, we can still span the codeword over more than one block to exploit diversity 
and achieve higher outage capacity if the power allocation and the codeword length K are smartly chosen. 




Number of blocks 



Fig. 4. Outage capacity v.s. number of blocks K, 1/A = 4, average power P = 3. 
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V. Outage Probability Over Parallel Dying Channels 

In the dying channel example of cognitive radio networks, secondary users have access to vacant 
frequency bands that are licensed to primary users. Some primary users may suddenly show up and take 
over some frequency bands, which results in connection losses if these frequency bands are being used by 
certain secondary users. Hence, each sub-channel (a frequency band) may have a different random delay 
constraint for information transmission due to the uncertainty of non-uniform primary user occupancy 
patterns. Specifically, the above system can be modeled as follows. Given a link with N parallel sub- 
channels as shown in Fig. [51 the codeword is spanned in time domain over K blocks and also across 
all the sub-channels. In some sub-channels, random attacks terminate the transmission before it is 
completed such that less than K blocks are delivered. For other sub-channels, K blocks are assumed 
to be safely transmitted. What is the maximum rate for reliable communication over such a link? For 
the single channel case, it turns out that there is no way to achieve arbitrarily small outage with a finite 
transmit power. However, in this section we show that an arbitrarily small outage probability is achievable 
by exploiting the inherent multi-channel diversity. 



T 



K 




1 2 3 N 



Fig. 5. Parallel dying channels. 



In this section, we extend the results of the single dying channel to the parallel multi-channel case. 
Definition 2: The outage probability of the parallel multi-channel case is given as 

f N L, ^ 

PoutiR, P^N) = Pi I Y,-Y1 + "?^/^) < ^ f ' (26) 
t i=i k=i J 

where R is the total rate over N sub-channels, a^.*"* is the fading gain of block k at sub-channel i, N is the 
number of sub-channels, Li = mm{K, [TJ } is the random number of surviving blocks at sub-channel 
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i, K is the number of blocks over which a codeword is spanned in the time domain, and P is the total 
average power such that P/N is the average power for each sub-channel. Since the asymptotic behavior 
is concerned, uniform power allocation is assumed over N sub-channels. According to different attack 
models, in the next two sections we investigate the asymptotic behavior of the above outage probability 
in two cases: the independent random attack case and the m-dependent random attack case. 

VI. Independent Random Attack Case 
Let the average power P be finite. Since log(l + x) ~ x if |x| ^ 1, when is large, we rewrite (|26l ) 



as 



t i=i k=i ) 



We assume that the fading gains s are i.i.d., and let the random variable Fj be 



k=l 



ii) 



For the case of independent random attack, we assume that L/s are i.i.d., and hence Y^'s are i.i.d. 
The outage probability given by (1271 ) can be recast as: 



Pout 



{R, P, A^) ^ Pr 1 1 ^ < i?/p| . (28) 



Since F/s are i.i.d., according to the central limit theorem, as the number of sub-channels A^ ^ oo, 
we have 

_Y,Y.^Af{fiy,a'y/N). (29) 

1=1 

According to Theorem 7.4 in [fT6ll on the sum of a random number of random variables, we derive the 
following relations: 

fiy = ^E{L)E{a) (30) 

4 = -^[E{L)Var{a) + Var{L)E{a)% (31) 
A^ 

where a is a nominal random variable denoting the fading gain, L is a nominal integer random variable 
denoting the number of surviving blocks of each sub-channel, and E{-) and Var{-) denote the expectation 
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and variance, respectively. As such, the outage probability can be approximated as: 

pUR.P.N)^^ ^'^~/^ ). (32) 

(Jy/V A/ 

As ^ oo, Yl!i=i converges to /iy. The outage probability decreases to over N if R/P is less 
than /iy, or converges to 1 if R/P is larger than /iyQ That is, even though all sub-channels are subject to 
fatal attacks, the outage probability can still be made arbitrarily small when A^ is large enough if the rate 
per unit cost is set in a conservative fashion, where fxy is a key threshold. This is remarkably different 
from the single dying channel case in which the outage probability is always finite since there are only 
a finite and random number of blocks to span a codeword. 

VII. m-DEPENDENT RANDOM ATTACK CASE 

In the previous section, we discussed the case where L^'s are independent. However, in a practical 
system, such as cognitive radio networks, the primary users usually occupy a bunch of adjacent sub- 
channels instead of picking up sub-channels independently. Thus, the Li's, across adjacent sub-channels are 
possibly correlated; and consequently the achievable rates across adjacent sub-channels are also correlated. 
On the other hand, if two sub-channels are far away from each other, it is reasonable to treat them as 

ll wit 



independent. Thus, we assume that Fj's are strictly stationary 
and variance. 



and m-dependentj with the same mean 



A. Central limit theorem for m-dependent random variables 

We first cite the central limit theorem for stationary and m-dependent summands from IfTSl (Theorem 
9.1 therein). 

Theorems (Hoejf ding and Robbins): Suppose {Xn,n> 1} is a strictly stationary m-dependent se- 
quence with E{Xi) = /i and Var(Xi) = < oo. Then as A^ — oo, we have 

1 ^ 

-= 5^(X, - /i) ^ Ar(0, vj, (33) 
* 1=1 

where t;™ = cr^ + 2 Y^Li Cov(Xt, X^+j) with Co\{Xt, X^+j) the covariance of Xt and X^+j. 

Proof: The detailed proof can be found in [[T9l . ■ 

^R/P is interpreted as the rate per unit cost in [17| . It is interesting to see that the quantity of rate per unit cost plays an important role 
here, which is due to the fact that we operate over both a finite power and a finite coding length. 

^Call a sequence {X„,n > 1} strictly stationary if, for every k, the joint distribution of (Xn+i, ■ ■ ■ is independent of n. 

'Call a sequence {Xn,n > 1} m-dependent if for any integer t, the cr-fields a{Xj,j < t) and a{Xj,j > t + m + 1) are independent. 
Simply put, Xi and Xj are independent if \i — j\ > m. 
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B. Asymptotic outage probability 

As assumed, the random sequence {Fi, F2, ■ ■ ■ , Yn} is stationary and m-dependent, and y^'s have the 
same mean and variance. Then the covariance is given as: 



\h\> m 

7(/i) — /iy \h\ < m, 

where /ly is the expectation of Kj given in (|30l ) and 'j{h) = E{YiYi^h). Meanwhile, 



(34) 



Vr. 



(35) 



h=l 



Due to the fact that the fading gains ■* and aq'^^'' are independent if p 7^ g or /i 7^ 0, we could easily 
obtain 'y(h) for \h\ < m, h ^ 0, as: 



1 
1 



E 



a 



(i+h) 



,p=l q=l 



E{L,L 



iJ-ii+h) 



(36) 



Assume that L,; and Lj have the same correlation coefficient p \f \i — j\ < m and « 7^ j. The correlation 
matrix is given as 



C 



1 p ■■■ 
pi p ■■ 

: p 1 p 




■■ 







p 



■•■ 

■ ■ 





p 



■• p 

p 1 1 



Note that the following main results can be also derived for other correlation matrices. 
Then (|36l) is simplified as 

l{h) = ^,{pal + pl\ 



(37) 



where p is a non-negative correlation coefficient, pi and cr^ are the mean and variance of the random 
variable L, respectively. 
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Substituting dSO]), dlB, and ([37]) into ([35]), we have 



According to Theorem |3l we have 



vm = 4 + 2m^^ (38) 
^'"'+#(l + 2r«,p). (39) 



N 



By simple manipulation, we have 



. V,n 

N . . 
1=1 



1 ^ 

-J2Y,-^Af{^Y,Vm/N), (40) 

where //y is given in (l30l) and is given in (l39l) . Hence, the outage probability for the m-dependent 
random attack case can be approximated as follows when is large, 

Pout{R,P,N)^^ ^^^lZJ^ ). (41) 
Vv„JN 

As we see from (l38l) that Vm > o-y, comparing (l32l) and (|4TI) . we conclude that the outage probability of 
the independent attack case is smaller than that of the m-dependent case given the same setting when the 
rate per unit cost R/P is less than fiy and the number of sub-channels is large. 

VIII. Outage Exponent 

As we learn from the previous sections, the outage probability over parallel multiple channels goes to 
zero as increases if R/ P < fiy for both of the two attack cases. In this section, we investigate how fast 
the outage probability decreases as N increases for both cases, which is measured by the outage exponent 
GOll defined as 

S{t) = lim -l"gPo-^(^>^>^)^ (42) 

N—*oo N 

where t = R/P. 

A. Independent Attack Case 

According to the results in [|20l . we could derive the outage exponent for the independent attack case 



as 



= sup{st- A(s)}, (43) 

s<0 
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for yt <to, where to = and 

A(s) := logE[exp(sr,)] 

= log My (s), (44) 

with My(s) the moment generating function of Yi. According to Theorem 7.5 in [[T6ll . we have My(s) = 
h{f{s/K)) where /^(z) and /(s) are the probability generating function of the discrete random variable 
Li and the moment generating function of the continuous random variable of* , respectively, t 

Example: If Rayleigh fading is assumed, af^ is exponentially distributed; hence the corresponding 
moment generating function is /(s) = (1 — s/Aq,)^^, where is the parameter for the distribution of 
the a), . Assuming that the random attack time has an exponential distribution, L is an integer random 
variable with following distribution: 

wq = Fi{L = 0} = Pr{0 < T < 1}, wi = Fi{L = 1} = Pr{l < T < 2}, ■ ■ ■ , 
wk-1 = Pr{L = K - 1} = Pr{/i - 1 < T < i^}, wk = Pr{L = K} = Fi{K < T}. 

Thus, we have h{z) = J2iLo'^i^^ ^^'^ My{s) = J2f=o'^ii^ ~ s/{K^a))~'^- Then we can derive the outage 
exponent numerically by solving (l43l) for a given t. 



B. m-dependent Attack Case 

For the m-dependent attack case, the techniques used in deriving the outage exponent for the independent 
attack case does not apply any more since here 1^'s are not independent. In this case, since the outage 

''The moment generating function of the sum of a random number of random variables, i.e., Sl ~ Xi + X2 + • • • + Xl, is the compound 
function h{f{s)), where L is a random integer independent of Xi, h{z) is the probability generating function of L, and /(s) is the moment 
generating functions of Xi. 
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probability has an approximate normal distribution, we have 





R/P-flY \ 



liY-R/P \ 



< exp 



2 



V / 



= exp I — 




) 



(45) 



Therefore, an approximate outage exponent can be quantified from the upper bound as 



if^Y - R/Pf 



(46) 



where is given in (l39l) . The outage exponent obtained by (l43l) is derived by using the large deviation 
techniques. Thus, it is exact while the outage exponent given by (|46l ) for the m-dependent attack case is 
approximate. However, when R/P <^ jjy, this approximation is accurate since the exponential bound is 
tight for the Q-function when its argument is large. 

Numerical results are provided here to validate our analysis for the parallel multi-channel case. We 
choose the random attack time T to be exponentially distributed with parameter 1/A = 5 and K is chosen 
to be 5. Rayleigh fading is assumed and the fading gain a^*^ is exponentially distributed with parameter 
1 and the noise has unit power. First, we demonstrate the convergence of the outage probability for the 
independent attack case and the m-dependent attack case, where the value of fiy according to the above 
simulation setup is 0.571. For the independent attack case, as shown in Fig.|6l the solid and dashed curves 
are derived by (|32|) while the circles and crosses are obtained by simulations. We also observe similar 
convergence for the m-dependent attack case in Fig. |7J In both figures, the outage probability goes to 
if R/ P < fiy, or goes to 1 if _R/P > /iy. We see that the accuracy of Gaussian approximations is 
acceptable with reasonably large values. 

Second, we compare the outage probability performance between the independent case and the in- 
dependent case. Here P = 2 and R = 0.5 nats/s. As shown in Fig. [H the outage performance of the 
m-dependent case is worse than that of the independent case even when m = 1 and p = 0.8. This is due to 
the fact that when R/P < /iy, the independent attack case is expected to have a smaller outage probability 
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X X-r^C X . 



analytical outage probability R-0.5 nats/s 

O simulated outage probability R-0.5 nats/s 

analytical outage probability R-2 nats/s 

X simulated outage probability R-2 nats/s 




O O O O o o o o o o 6 

15 20 25 30 35 40 

Number of sub-channels N 



Fig. 6. Outage probability convergence behavior of the independent case: /iy =0.571 and P=2. 



X I X X X jj ! ,X'-X--X- ? < " X X-i - X - X -X K X i X 



analytical outage probability R-0.5 nats/s 

O simulated outage probability R-0.5 nats/s 

analytical outage probability R-2 nats/s 

X simulated outage probability R-2 nats/s 
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Number of sub-channels N 
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35 40 



Fig. 7. Outage probability convergence behavior of the m-dependent case: /i-i'=0.571, m=l, p=0.8, and P=2. 



as we discussed at the end of Section IVni-B[ However, the outage probability of the m-dependent case 
still decreases to but at a slower rate as the number of sub-channels increases, which is caused by 
the fact that the m-dependent attack case has a smaller outage exponent. 

In Fig. m we compare the various outage exponent values between these two cases over the rate per 
unit cost R/ P with the simulation setup as follows: A' = 5, m = 1, and p = 0.8. First, we see that 
the outage exponent for the independent attack case is larger than that of the m-dependent attack case 
when the average attack time 1/A is the same. Second, for both of the independent attack case and the 
m-dependent attack case, a larger average attack time 1/A results in a larger outage exponent. 
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O independent attack case 
V m-dependent attack case 



o y 'V 



o o 



o o 



10 15 20 25 30 35 

Number of sub-channels N 



Fig. 8. Outage probabilities comparison. P=2, R=0.5 nats/s, m=l, and p=O.S. 




0.1 0.2 0.3 0.4 0.5 0.6 07 

r=R/P 



Fig. 9. Outage exponents for independent and m-dependent random attack cases: m=l, p=0.8, and K=5. 



IX. Conclusion 

In this paper, we considered a new type of channels called dying channels, where a random attack 
may happen during the transmission. We first investigated a single dying channel by modeling it as a 

-block BF-AWGN channel with a random delay constraint. We obtained the optimal coding length K 
that minimizes the outage probability when uniform power allocation was assumed. Next, we investigated 
the general properties of the optimal power allocation for a given K. For some special cases, we cast 
the optimization problem into convex ones which can be efficiently solved. As an extension of the single 
dying channel result, we investigated the case of parallel dying channels and studied the asymptotic outage 
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behavior with two different attack models: the independent-attack case and the m-dependent-attack case. 
It has been shown that the outage probability diminishes to zero for both cases as the number of sub- 
channels increases if the target rate per unit cost is less than a given threshold. Moreover, the outage 
exponents for both cases were studied to reveal how fast the outage probability improves over the number 
of sub-channels. 

APPENDIX 

A. Proof of Theorem [7] 

Let us consider minimizing the outage probability given by (fTTl) . When K = 1, the proof is trivial. 
When K = 2, the outage probability is 

Pout 

(2) = Wo + Pr{log(l + aiPi) < 2R}wi 

+ Pr{log(l + aiPi) + log(l + a2P2) < 2R}w;. 

As we see from the above equation, if Pi < P2, we have Pr{log(l + a;iP2) < 2-R} < Pr{log(l + aiPi) < 
2R}. Hence, we can achieve a smaller PoutC^) by swapping Pi and P2, since the last term in PoutC^) is 
not affected by such a swapping while the second term is decreased. 

When K > 3, for any j > i, G {1, ■ ■ ■ , K}), if Pi < Pj, by swapping Pi and Pj, all the terms 
containing both Pi and Pj, i.e., all the probability terms in the form of Pr{- ■ ■ + log(l + ajPj) + ■ ■ ■ + 
log(l + ajPj) + ■ ■ ■ < KR} will not be affected. However, the probability terms containing P/ but not Pj 
can be decreased by such a swapping. Thus, we could achieve a smaller outage probability in total. 

Therefore, the optimal power allocation profile over i.i.d. fading is always non-increasing, i.e.. Pi > 
P2 > ■ ■ • > Pk > 0. 

B. Proof of Theorem |2] 

When the coding length K = 1, the outage probability is 

Poutil) = Pr{log(l + aP) < R} Pr{T > 1} + Wo. 
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When we choose any other arbitrary values for K, i.e., K = M and M ^ 1, according to (fTTl ). the outage 
probability is 

Pout{M) 



+w;o + 5^ Pr i ^ ^ log(l + aPi) < r\ w,. (47) 

i=l I 1=1 J 

Due to the concavity of the log function, we have jj J^fii log(l + (^Pi) ^ log(l + ctP). Hence, 

Pr I ^ ^ log(l + aPi) < i?| > Pr{log(l + aP) < R}. (48) 

Moreover, it is obvious that summing over only a portion of the M blocks yields an even smaller value, 
i.e., jj log(l + ^Pi) ^ log(l + "^)' with 1 < i < M - 1. If 3Pj > 0, for i < j < M, the strong 
inequality holds. Therefore, we have 

Pr I ^ ^ log(l + aPi) < i?| > Pr{log(l + aP) < R}. (49) 

Noting that J2iLi^ '^i — Pr{l < T < M}, and considering (l48l) and (|49l) . the following inequality can 
be derived for (l47l) : 

Pout{M) > wo + Pr{log(l + aP) < R} ■ 

(Pr{T > M} + Pr{l < T < M}) 
= Pout{l). (50) 

From (|50l ). we see that Pouti^) has the smallest outage probability when fading gains are the same, which 
means that the optimal coding length is K = 1 with Pi = P. 

C. Convexity of the optimization problem in f l27]) 

We first check the Hessian matrix of the objective function in terms of Pi. 



The jth term is: 



(51) 
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^ , \ 

1 

p2p2r-rJ-l p. 
-^2-^j lli=l,i7^2 

Let z G M^^, then 

where p(^) = (1/Pi, I/P2, ■ ■ ■ , 1/^,, 0, ■ ■ ■ , O)"", and M = diag (^^, ^, ■ ■ ■ , 0, ■ ■ ■ , o) . 

Therefore, (|5T1) as the summation of all the A' terms is positive semi-definite. Hence Pout is a convex 
function in terms of Pk- In addition, Pk lies in a convex cone as shown in Theorem. [U Hence the 
problem is a convex problem. 

D. Sufficient conditions for the convexity of the optimization problem in d25l) 
Let /i : — >W., g: M" — ^ M\ and / = /i o ^ : M" — ^ M be defined as: 

/(x) = h{g{x)), dom / = {x G dom G dom h}, 

where dom is the domain of a function. The following two lemmas can be established. 
Lemma 1: / is convex if h is convex and nondecreasing, and g is convex. 

Proof: See Chapter 3 in IH. ■ 
Lemma 2: The outage probability function given by (|24)) is convex, if 

KR- log Pi < 0, 
ATP -log Pi -log P2 < 0, 



2 



K 

ATP log Pi < 0. 
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Proof: As we know, the error function erf(x) = e * (it is a convex function for a; < and it 

is non-decreasing. 



Since gmiPx) is convex and erf(x) is convex for x < and nondecreasing, according to Lemma 1, if 
gm{PK) < 0, erf((7m) is convex over Pk- Hence, the objective function given by (|24l) is convex under 
the conditions given in the lemma. 

Since the constraints in (l25l) are obviously convex and the objective function is convex under these 
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